[ 149 ] 



V. A New Treatment of Optical Aberrations, 



By K. A. Sampson, F.R.S. 



Eeceived March 15,— Kead May 23, 191 2^. 



The method developed by Gauss in his * Dioptrische UntersuchTuigeii/ is probably 
the most powerful, as well as the readiest, method in geometrical optics. It has 
in effect hitherto been restricted to systems in which the relations of original and 
emergent rays are strictly linear, or, in optical language, those in which the 
aberrations can be neglected. It is true that Seidel bases his celebrated discussion 
of aberrations upon Gauss's method, but he soon modifies it and replaces its system 
of co-ordinates and characteristic steps by others. The following pages show how the 
method may be extended and retained throughout the discussion of the aberrations of 
any co-axial system. They will be found to throw light upon the general relation- 
ships of the well-known Petzval condition and Abbe Sine condition, to furnish a ready 
method of describing, analysing and measuring the faults of an optical image, and to 
be particularly adapted to numerical calculations, to the order to which these are 
necessary for telescopic objectives. 

It will be convenient to state here the essentials of the method in the form in 
which they will be used later. Let Oxyz, OV^/V be rectangular axes in the original 
and emergent media, of which the refractive indices are ju, jJ respectively. Ox, O^x 
are the axes of the optical system. Take the equations of any ray before and after 
its passage through the system in the respective forms 

y = fix-j-h, z = yx^-c, 

and ........( 1 ) 

then, provided there is a strict linear correspondence as well as symmetry about the 

axis, we may put 

V = gb-^ h^, d = ge 4- Ay, 
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where </, h, k^ I are constants involving the curvatures of the refracting surfaces, the 
distances between them and the refractive indices ; also 

gl — hk = /ii//ii\ 

FoUow^ing Seidel, w^e shall call such systems normal systems. 
In particular, for a single refracting surface, 

2x — B(?/ + 2^^)-f...5 



without change of origin^ the scheme 
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where ^ is put in place of zero. Or again, a simple shift of origin by a distance d 

may be represented by the scheme 
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If two instruments be represented by the schemes 
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light passing through (l) first and then through (2), and the emergent origin for the 
first being made the same as the original origin of the second, their combined effect is 
given by the scheme 

"9192 + Kth^ hg2 + Kh 2 " 



< 



>. 



gjc,^ + k^k, Jhh + hk 



(3) 



-^ 



which may be written down by multiplying the rows of the later scheme into the 
columns of the former, as if they w^ere determinants. It will be shown hereafter that 
this rule is remarkably well adapted for numerical calculation — a fact that does not 
seem to have been remarked before. The scheme corresponding to any system, as, for 
example, any thick lenses, arranged at intervals along an axis, may be built up from 
its elements by this rule, by writing down the schemes in order belonging to the 
successive refracting surfaces and shifts of origin, and compounding these ; if we have 
to coinpound in this manner a sequence of schemes 
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lien, provided we do not change the order in which the schemes present themselves, 
the composition may be effected in such groups as may be convenient, and may be 
performed either from left to right or from right to left.* 

The product of the determinants of the component scheme (gili—hiki)^ {Q^h'^^h^^)'^ ••• 
gives the value of the determinant GL— HK of the compounded scheme. 

The analytical scheme corresponding to any instrument of which the cardinal 
points are known may be w^ritten down at sight, and conversely, by the relations 



»■ ■ . ■ .« ■ «■■"« ■ '" ■ 

O F 



H N 



O' H' H 
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Fig. 1, 



J/T?/ 



iW^ljh HF = r#, NF=:1/^; O^W ^ ^glh HT^-^l/1, Wl 



n/k, 



where 



n = gl—hlc . . . . 



(4) 



and H, H^, N, N^ F, W denote, as usual, the unit points, nodal points, and principal 
foci respectively. 

If it is desired to work geometrically, we may set the original and emergent axes 
across one another in a figure at any angle, H and H' being superposed (or else 



f\. 



xi m li iwHn \t\tm 



H 



P\ 



f 



H' 



N 



Fig. 2. 



W), and finding a point with co-ordinates HF, WW ; then any straight line 
OPP' through determines points PP^, which are conjugate foci. 

The following method of compounding any two given systems may also be 
mentioned : — - 



* A general discussion of the linear system by the author will be found in 'Proe. London Math. Soc.,* 
vol. 29, p. 33. 
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Set the axes as shown in the figure, the distance between the points HiH/j and 
HgH'g being annulled. Then the lines FgQi, F^Q^ give #, #', the principal foci,, and 
the line OiO^ gives points N, N^ which are conjugate to one another and are the nodal 
points of the compound system. 

We see that it is always possible to determine a geometrical system that shall 
correspond to any given values of g, h, k, L Thus, for example, n == implies that 
F' is conjugate to evef^y point of the original system, or, what is the same thing, that 
every emergent ray goes through F'. 

If the emergent origin is at the principal focus, g ™ 0. 

If the original origin is at the principal focus, ^ = 0. 

If the original and emergent origins are conjugate points, h = 0, 

We shall now consider the case of refraction of a general ray at a symmetrical 
surface centred upon the x-Sixis and shall show that a scheme {g + Sg, ...} maybe 
derived for it, which shall include the aberrations ; these, represented by the 
additional terms Sg, ..., will, of course, vary from point to point with the squares and 
products of the co-ordinates and angles of incidence upon the surface, whereas for the 
pure linear scheme g, ... are the same for every ray of the beam. 

Taking rectangular axes Oxyz, let the equation of the surface separating the region 
of index ^ from that of index ij! be 

2x = 'B{y'-^z')-^\Q{y' + zJ, ,....,., (5) 
Let a ray 

y =^ I3x + h^ z = yx + c, 

in the original medium be transformed by refraction at the surface into 

y^ =■ ^^x^ -f ?/, / = yX^ -f- c\ 

where the axes are, in fact, the same but are accented to indicate the difference of 
medium. 

The positive direction of the .^-axis is that in which the light is travelling. 
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In the diagram Pq is the point where the original and emergent rays meet at the 
surface, PPq is the original ray, PTq the emergent ray, and (O, &, c) are the co- 
ordinates of P, (O, b\ c') those of P^ and we shall take (^o? ^o? ^o) ^s those of Pq. 




If {I, m, n) {l\ m\ n') are the direction cosines of an original and emergent ray, 
(^, g, r) those of the normal to the surface at the point of incidence, we have the 
known equations 

{ijI'-'iJV)Ip = {iixm—ij!w!)lq ^ {fj^n—ijln^^lr = jm cos O—jul^ cos 6\ 

where 0, 6^ are the angles made by the two rays and the normal. 

Now 



l = m/^ -n/y -l-i^'-h' 




I'-m'l^' -n'ly' -1 IP" ly" 




p = g/B6o+iC6o(V+^«^) = r/Bco+iCco(&,/+Co^) = -l+iBV+iBV = 


-l+|-gH|r 



if we neglect higher powers of the small quantities. 
Further 

cos = 1 -i0^ cos e' = 1 ~|0'^ 

where 

and we have approximately 



-tY 



M (^-g) = / {^'-ql /x (y-r) = m' (/ ~r). 



Substituting above for m, m\ we have 



^P{l-iP'-h')-f^'^'{l-W'-h") = {f^-f^')q-iq{f^e'-f.'e'% 



But since 
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6/ JUL = 0//X. 
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therefore 
and 



ee' = (^_g)(/3'-g) + (y-r)(y-r). 



Hence the right-hand member of this equation i^eads 

+i (m-mO (A5'+yyO-i (i8+/3') {Mi8-M'/30-4 (y +y) (mt 

= g[(M-/)+i(M-M0(g^+^^)-iM(/3Hy^)+V(^'^+ym 



-mV)] 



or, since 
and 

therefore 






and the equation becomes 

or dividing by the coefficients of jS' and writing 



a 



/m^ 



rL 



^^ = 6[^(l^,^)B~-|(l-tOC(6^+c^)+itiBO^^ + y'-'^'~-y')] 



Also 

Therefore 
but approximately 
therefore 

b 

or if we write 



we may put 



+ /3[n+in(/3'^+y'^-^^~/)-i(l-n)(gHr^)]. 

bg = h + Uofi = h' + a„^'. 

h' = b + a,{^-l3'); 
13' = -(l-n)B6+«/3; 
^b[l+^{l-n){q^ + r^)]+/3[i{l-n){g^+r^)/B-]; 

a, = ^{l-n){q' + r') = i{l-n)B'{¥ + c'), 
b' = 6[l+w] + /3[w/B], 



» • s s » s 



(6) 



^' = 6 



— (l— n.) B + B\jr— t^ ft) 



+/3[ra+^— ft)]. 



t • • 



(7) 
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In the same way it follows 

c' = c [l +0)] +7 [ft)/.B], 



y = o 



±5" 






for the case of the paraboloid 

C/B^ - 0, 
sphere 

C/B^ = B. 
We shall generally write 

C/B^ = 6B. 

We remark that the coefficients that transform the (&, /3) system into {b\ /3^) are 
the same as those which transform (c, y) into {g\ y), and for any surface each is 
expressed in terms of the two functions ^, w defined by equation (6), in addition to 
the refractive index and curvature. 

These equations therefore permit us to treat rays which cross the axis with the 
same readiness as those which intersect it, a thing which is very troublesome in the 
trigonometrical discussion of the question. They also apply equally easily to the 
sphere, the paraboloid, and any intermediate form. 

Before proceeding with the discussion of these formulae I shall verify that they 
cover the known expression for longitudinal aberration on the axis after refraction at 
a single spherical surface, as it is given in the text books. 

Suppose the ray meets the axis Sit x = v^ x^ ^ v\ so that 

h-\-^v = 0, 6'+/3V = 0; 
then the equation connecting v^ v' for the case of the sphere is 

— W[-(l--n)B + B(^/^--w)]+'?;^[ti+^--ft)]---i^[l+w] +ft)/B = 0, 
or, dividing by vv' and rearranging the terms, 



but 



^-^-(l-^^)B:=(B-i)(-V^H-.^:^); 
V V \ v/\ i5v / 
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Hence the right-hand member above is equal to 
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This is one of the usual expressions; compare Herman's * Optics/ p. 189, (iii.)- 
After a slight transformation it leads to the Zinken-Sommer expression for the 
separation of the focal lines in any co-axial system, and thence, as Whittaker has 
shown (* Theory of Optical Instruments,' p. 26), to the expressions of Seidel's 
theory. 

We may verify also that these expressions lead to the known results in the case of 
the parabolic mirror. 

Consider the focus for rays parallel to the axis, ^*.e., when ^ = y = 0. 



But 



so that 



== i+a)+'y'[-{l-r6)B + B>/r]. 
v' = ll-^B'{¥+c')]/[2B + 2B'{b' + c')l= 1/2B, 



so that the longitudinal aberration vanishes at the principal focus. More generally, 
the ray y^ = ^'x-^-V] ^ = y'x-\-c\ corresponding to a general incident ray for which, 
say, 7 = 0, meets the plane x^ = d^ in points whose co-ordinates are 

y^ =: &[-~.(l-9^)Bc^'4•l+B(i''^+ft)]+iQ[(i'(n+\^--w)+fe)/B], 
z' = c[-{l-n)Bd'-^l+Bd'yfr + (ol 

For the paraboloidal reflector 
since /3' = — ^8— 2B6, y = 2Bc ; therefore taking the focal plane d^ = l/2B, we have 



y^ = -/3/2B~B6^^-/3[|-B(fe^ + c^)+6/3], 



z = 



-B6c/3 



...♦•. \*^ / 



these are known expressions, leading to the theory of the coma of a parabolic 
reflector; cf, Plummer, ' Mon. Not. E. A. S.,' LXIL, p. 365, (9), (10). 
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If we compare the scheme 
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— (l— '^)B+B^— eBw, tn-^ 



ft) 



. • (7). (8) 



with the equations (4), we see the aberrational scheme is the equivalent of a linear 
scheme, for a single surface, for which the plane of origin passes through the point 
X = m/{l—n)3^ that is, through the actual point of incidence £c = |-B(6o^+Cq^)...; 
if we take the curvature as B (l +|-€g^+|-er^), and the refractive indices ^ (1— i-i8^--iy^)? 
/mf {l—^^^^—^y^) respectively. 

As to the latter it may be noticed that the exact equation for refraction of a ray 
impinging at the origin, and in the plane Oxy^ is 



or, to our order. 



jii sin {tan~^ ^) == /i' sin (tan^-^ j8^) 



Hence the aberration i/r may be described as due purely to the obliquity of the ray 
to the axis, and the aberration w to the lateral separation from the axis, and we see 
that the somewhat remarkable fact that two functions ft?, %fr suffice to express the 
aberrations of every ray may be stated in the form that there is no term which is 
produced jointly by obliquity and lateral separation. 

If in any instrument we have a number of surfaces each introducing aberrational 
terms, and if the schemes preceding and following the surface (r) be compounded so 
as to read, say, {g, ...}, {g\ ..,}, then the whole may be represented by 
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and the portions added to the general scheme in consequence of the aberrations of 
the r*^ surface will be 
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(11) 



and in this the schemes {g^, ...}, {^, ...} may be taken at their "normal"' values 
without regard to aberrations introduced by surfaces other than the surface (r). If 
then we write, adding the effects of all the surfaces, 






(12) 



and if we now denote by {^, .-.} the scheme got by compounding all the normal 
schemes of the instruments in succession, whether these are refractions, or mere shifts 
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of origin from one surface to the next, we have for the relation between any original 
ray 



and the corresponding emergent ray 



y^ = /B^x' + h^^ z^ = y V 4- e\ 



[k+sk, l+sl 



> = {b' + Sh', ^'+S/3'), 



(c, y) 



Gr + oGr, Jtl + oH. 



• « * \ JL %J f 



> ^ {c^ + Sc\ y' + dyO: 



where b\ ^\ c\ y are the values that would result if ^G, ... were all zero; and 
turning back to the expressions (6) we see that co,,, -i/r^. are expressed in terms of the 
incidence — (6, ^, c, y)--—upon the plane Oyz by the equations 



where 



(^r — 2 



I ^ ._ II; \ _Oj, \ C/,^, "t" '^J» /5 



where 






. (14) 



^^ = kj} + l^^ 



y,. = k^e + l^.y, 



^'r = K+i^+h^A Yr = h+i(^+K+ir 



We notice that the determinants-™(gZ - M)— of the scheme multiplied by each t/^^ 
will always be zero, and that of the scheme multiplied by w,. also, for the case of the 
sphere. This supplies a useful check. 

The numerical management of these formulae for actual systems is dealt with later. 
I shall now consider their analytical and geometrical properties. 

If 

fG+^G, n+m 

(&,/3)< 

I^K + ^K, L + ^L 



> = {h', n 



and SG, ... are quadratic functions of b, c, /3, y, with the symmetries imphed in the 
forms above, we may put 



m = il^iH {b'+c') +2S,K (b^+Cy) +4H {/3'+y')}, 
SK = i {S,K {b' + c') + 24K (6^ + cy) + 4K (/3^ + y')}, 
SL = i{8Ju{b'+c')+2S,L{b^+cy)+§,L{^'+y')}. 



■ (15) 
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In these expressions the values of ^iG, ... are not unrestricted. Thus, for example, 
the rays which originate in the point (6, c) must upon emergence be normal to a 
surface. Consider the conditions that 

where V , /3', c\ y are functions of two variables ^, y, as above, should be normal to a 
surface. 

If {x\ y\ z^) is a point upon the surface, then we have for all directions upon the 
surface 

ddd -\- dy^ -\- y dz^ = 0, 

or, since x\ y\ %^ are functions of j8, y only, 



Also 



Therefore 



a^ ^^ ^/3 ^^ 3^ ■" ^' 3y ^^ 3y ^^ ay "~ ^• 

yf = ^V + 6'5 ^ = yV + c^ 

a^ 

By \ 3y 3y/ \ ay ay 



(l+^'2 + y2) + ^/L/M + p;\^yL/Sz + ^) = 0; 



or, say. 



A {a:' (1 +^VyTn + [f I' +y ^) (1 +/3'^ + y'^)-'^= 0, 
^ (a;' (1 +^'^ + y^)i/^} + f^' 1^ +y 1^) (1 +;3'^+y2)-iM = 0, 

dy \ dy dy / 

so that the necessary and sufficient condition is 

3 (Ay) 3 (Ay) 

Retaining only the terms of lowest order we have 



^ij 



ll = H ^ = 4G . 6c + 4G . 6y+^2H . c/3+<53H . /3y, 

dp dy 

TTTT = L ::^ = 4K . 5c + 4K . 6y + 4L . CiS+c^oL . By, 

op dy I o I /, 

II = 4G . hc^-lf^ . c^+4H . 6y+4H . /3y 1^ = H, 
^ = ^K . &C+4K . C;8+4L . &y+,53L . /3y ^' = L, 

cp dy 
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wlienoe the condition 



H 4Gr— 4H I (6y— C/3) = 0, 



Ju OoJA. — Ooi-J 



OoVJT — OoJOL OoJX. — OoIj 



'3 



3-^ 
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^^?, say. 



» • • .• » . . ♦ I X VJF f 



Also, this result must remain valid if we pass the emergent beam through any 
further optical system. This is a step that must frequently be taken, and it will be 
convenient to write down generally the formula to which it gives rise. 

If we have 



< 



g+Sg, h+Shl fg'+§g\ h^ + Sh^^ 



< 



K-hSK L + ^L 



> 



• n 



and if 



Sg = i{Sig{b'+(^)+2i0{^^^^ ..., 



J > 



and, further, 



SQ = |{(!iG{6^+c^) + 24G(&/3+cy) + 4GO^+y^) 



i ?. • • M 



then the following formulsB result :- — - 

= g%g + h%k + g {gV + 2g^4^' + m^' ) + k {g\U + 2gU,h' + h%h'] 



8iQr - 

\'E:=g'SJi+h'Sjl+h{ 



ibid. 



c -^ I 



J 



ibid. 



}, 



k%g + VSM9{g%k'+2gkd,¥+ mj(/}-¥k {g%'¥ + 2gkS,F +1^4^; 
MSih + V Sil+h{ ibid. } + 1 { ibid. } , 



4^ 

U2XX 

Oo-*-^ 

03IJ 

OgXX 



g%g-^h%k+g {ghhg'^igUhk) §29'+klS,g^}+k{ghAh'+{gl+hk) SJi^^MSJif], 
g^§2h + h^S2i +^{ 'ibid,:- } i-l { - - ibid. }, 

k%g+V S,k+9 {ghS,¥+{gUhk) S,y+kU,M} +k {#J/ HgUhk) S,V -^kU,V }, 
MS2h + VS2l+h{ ibid, }+M ibid. }, 

g^s9'^h%k+g{¥.S^%2MS29' + l%9'}+k {m^h^ + 2MS2h^+m^'}, 
g^SJi-^h'S^l +h { ibid. }+? { ibid. }, 

F4^+l'4^+/M, ^*&^*<^. }+M ^*^*'^- }• 



These formulae with 



G = gg^-hkh^. 
K=^gM-^M\ 



== fig -r (jth 5 
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are of fundamental importance and cover all cases; they will be quoted as 



> • • • 



(17) 



Apply them to the equation (16); we have 

4G-4H = g' {S^-S,h) + h' {S,k-S,l)+hn {S,g'-m') + lnW-S,h'), 

S,K-S^ = k'{S^-S,h) + l' {S,k^S,l) + hn{S,k'-§,r) + ln{S,k'-S/), 
where 

n = gl—hk. 
If we write in this 






we have 



S,G-AIl = p{hg'+lh')+n{h{S^'-SJi')+pW}, 

= {p+np')K+nh{{S^'-SJi')-)?'g'}, 
S,K-S,L = {p+np')L+nh{{§,k'-S,h')-p'k'}. 

In the same manner we find 

S,G-SJI = {p+np')G+nk{{S^'-S,h')-p'h'}, 
S,K-SJ. = {p+np')K+nk{{S,k'S,l')-p'l'}. 

Compare these with (16) and remember that the two systems (gh ...), {g'h' ...) 
are arbitrary and independent of one ainother. Then we see that if for these systems 

S^g-S^h S^-§,h _ S,k-S,l _ S,k-S,l _ . 
g h ~ k ~' I ~^' 

Kg'-Kh' _ _ _ _.,, 

g 

^^G—S^rL _ ^aG — ^2-H. _ ^2-^ — ^iL _ (^sK — ^2-L _ sn /i q\ 

Jtl xv ij 



then 



where 

Now If we examine the case of the single surface, for which 

8^ = (1 -^) B^ ^29^ = 0, (^3^ = 0, ^i/i = (1 ^n) B, ^,h = 0, 8Ji = 0, 

and 

gr = 1, h = 0, k = — (l-!?^)B, I = n,- 

the conditions are fulfilled and we have 

^ = -(i-».)b = m(V--)b; 

VOL. CCXII. — A. Y 
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for two surfaces Bqj B2 with refractive indices^ ijl__i, fXi^ /ul^, as in Sbidel's convenient 
notation 

P + np^ = M-i [(-i -— )Bo + (i ~^1) B, 
and for any sequence of surfaces whatever 



vA*2r + l ^2r-~V 






This will be recognized as the expression which figures in the well-known '' Petzval 
condition for flatness of field." It was given by Petjzyal without proof in 1843, and 
it is a comment upon the difficulty which the geometrical method finds in removing a 
condition that may have been tacitly introduced that its proper position has so far 
remained obscure. Its general geometrical implications will be considered later. 

Besides the condition that the rays of any thin bundle should always be normal to 
a surface there is another general property to which they are subject in all systems. 
For normal systems in which we have stigmatic correspondence this is usually called 
the Helmholtz magnification theorem connecting the linear and angular magnifi- 
cations. For aberrational systems it would at first appear as if both linear and 
angular magnifications lost their meaning, but I have succeeded in generalizing the 
theorem in the paper already referred to."^ In the first place focal lines in the 
original system are shown to correspond one to one and not pair to pair with focal 
lines in the emergent system ; and rays which issue from any point in a focal line in 
a plane perpendicular to that line lie in a plane in the emergent system perpendicular 
to the conjugate focal line which they meet in a point. Such planes are called planes 
of correspondence. The behaviour of any ray may be traced through the behaviour 
of its projections upon the planes of correspondence. The separation of two parallel 
focal lines compared with the separation of their two conjugates preserves the idea of 
linear magnification and the angles in the planes of correspondence that of angular 
magnification. Then if a is the separation of two focal lines which lie parallel to one 
another in a plane perpendicular to an original ray at any point and a^ that of their 
two conjugates, and if a is the angle between two rays issuing from one of these lines 
in a plane of correspondence perpendicular to both and a^ the angle between the 
same rays on emergence, it is proved that 



jmaa = iji!a^a\ 



This is completely general. Now return to the case of surfaces centred upon an 
axis. It is clear that for any point ofp the axis, say the point (O, h, O), one of the 
planes of correspondence, is the meridianal plane passing through the axis and the 
point itself. 

* ^Proe. L. M. 8.,^vol. 29, p. 70. 
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Now consider the substitution 

ft'-(G4-(5G)&+(H + ffl)/3, 
/3' = (K + (5K)64-(L+5L)/3. 

Then if we shift the origin in the emergent system to d! the first will read 

Choose d^ 80 as to make the coefficient of ^ zero ; then 

V - {(G+^G)--(K + ^K)(H + (5H)/(L + ^L)} 6, 
and the coefficient 

G+^G-(K + ^K)(H + ffl)/(L+^L) 

is the linear magnification for narrow pencils emerging in the general direction (/3) 
from the point (O, 6, O) in the meridianal plane. Again from 

the angular magnification for the same is L + ^L. 

Hence (G4-(5G)(L + ^L) — (H + ^H) (K + ^K) is equal to the ratio of the effective 
refractive indices. But we have seen on p. 157 that the change of ray effected by an 
aberrational system is equivalent to the use of refractive indices /i(l--|-^^), 
jjL {1—^0^)^ ... throughout. So that the expression above is equal to 

where may be taken as the final value of /3 after any number of transformations ; 

or equal to 

N{l+i(K6+L/3)^-i/:f}. 

Identifying term by term with the expression above we have the relations 

4N = G^^L+L^G-H^K-^K^H == KLN, (20) 

4N == G4L+L4G--H4K-K4L == (P--l)N. 

The relations (20) may also be proved from a sequence formula out of the equation 
(17); thus 

JN n V n n n } 

^-§ = ^-^ + \gh^-4^igWhk)^4+uH], (21) 

N n I n n n \ 

JN. n I n n n 

Y 2 
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in fact, it was by such a method that I found them ; but their real significance, is 
contained in the proof given above. 

We have thus found among the twelve aberrational coefficients six relations which 
may be expressed in terms only of the focal length and other cardinal elements of the 
normal system, or seven, if we include the Petzval expression as of that class. 

Let us consider next what geometrical description can be given of the occurrence 
or absence of the twelve coefficients. It must be remembered that for different 
choice of origins the coefficients do not preserve an identity. Thus if we shift O, the 
original origin to the point (— c?, o, o), the new set — -^iG, ... — -is given in terms of the 
old set — "SiQ^ ,,. — ^by writing in the equations of p. 160. 

§ig = S^h ==.., = 

and 

g = ij h = d^ k = 0^ 1=1; 

and if, on the other hand, we shift the emergent origin to d\ we have ^iG, ... 
connected with SiQ, . . . , which now figures as the old set, as if in the same equations 
we wrote 

S^g' = §^h' = . . e = 0, 

g^ ^ 1, // = d'^ y = 0, V = 1, 

and in the event of both these changes being made a system (g+^Sg, ...) is trans- 
formed into (G + ^G, ...), where 

G — g + dlcy H = A + dg + d^l + ddlh, 

K = y^, L = l + dk; 

(5iH = SJi + dS,g + d'{Sil + dSik), 

§,K = S^k, 

SJui = §J, + dSik, 

(%G = S^+ dS^g + cZ'c^aK , 

4H = S,h +d{§,h+ S,g) + d\g + c^'^L, 

S,lj = U+d{\l+^Ji)+d\k, 

4G = §^gr\-2d§,g + d%g + d'l^, 

4H = SJi+ d{2S,h+S^g)+d' {S,h + 28^) + d%g + cZ'^L, 

4K = ^Jc+2d\k + d\k, 

4L = 8,Ud{2U-^8Jc)-^d'{U-^2U)^d%k. 
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But let us defer discussion of these and examine two particular cases of special 
importance, namely, let us assign meanings to ^^G, ... : (l) where the emergent origin 
is the principal focus, so that G = 0, and therefore ^gGr = ^iH, and (2) where the 
original and emergent origins are conjugate, so that H = 0, and therefore SJoc = ^H. 
In the former the original origin may be anywhere, but may conveniently be supposed 
to lie at the tangent plane to the first refracting surface. The original rays are in 
constant direction, so that we may take 

^ = const, y 7 = 0, and, say, h = d cos ^, c = d sin ^. 






Then if we receive the emergent ray on the plane parallel to O^y^z^ which passes 
through a point slightly removed from O', say at 

and it cuts this plane at j/ = b' + Sb', z' = c' -^^c\ we have 

b'^W = [* + K^/'] &+ [H+ W/'] /3 

+U cos ^ [d?\Q+Ufi cos .^^sG+^^^G] +i^ [(^^<5iH + 2(^/3 cos <^^2H+/3^^3H], 

+^dsi.xi<l>\d%G + 2d^(ios<t>8^G+0^§0]. . ... . . . . . . (23) 

Let us take 

b' = (K+LSf) 13, c' = 0, 
so that 

Sb' = i^ld'{SJI+S,G)+^%K] 

+ coB,pd[K§f'+y%G+il3'{S,G+2S,lS.)] 

+ cos 2^c?^ [i;84G], ^ 

■ Sg' = sin <i>dlKSf'+id%G+i^%G] 

+ sin 2<t>d' li^S.G']. . . ... . . . . . . . (24) 

These express the amounts by which the aberrations disturb the ray from its 
normal focus. Consider the lines in turn and examine their significance when the 
original ray traces out a circle d = const. 

The terms 

i^[d'{S,R+S,G)+^%K} or i^[2d%G+fi%B.'} 
give a fixed point. It may be considered as adding to the focal length 

the terms 

d'S,G+i^%R, 
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of which the former may be called the comatic increase of focal length and the latter 
the distortional increase of focal length, since these are evidently their characters. 
The two terms 

cos^d[K(5/^+|-d!^(5iG+|-^^(4G + 2(5,H)], ^ sin ^c^ [K^/^ + ici^^^G + i/^^^G] 

represent an ellipse which may be varied by choosing ^f^ at different values. If we 

take 

K.5/*^ + U%G + W (4G + 24H) = 

the ellipse becomes the primary focal line ; if we take 

K^/'+i(^^^iG+i/8'4G = 

it becomes the secondary focal line, in advance of the primary line by the amount 
/3^4H/K. Generally I shall call it the/oca^ ellipse and, as a rule, shall take 

which gives the focal circle 

d cos ^ [iiS^^H], -d sin ^ [i/S^^H] 

situated midway between the focal lines. This circle is described backwards as the 

original circle d = const, is described forwards. 

Finally the terms 

d' cos 20 [i-^^G], d' sin 2^ [i/^^G^ 

give another circle which I shall call the comatic circle ; its radius =^(3x comatic 
increase of focal length, so that they vanish together. As the original circle 
d = const, is described once, forward, it is described twice, forward, each point upon 
it corresponding to two diametrically opposite points of the original circle. 
Consider the focal circle and the comatic circle simultaneously ; we may take 

I GOB (/)-{- m COB 2(p, — I Bin (j> + m mi 2(1), 

where 

l = ^d^%Il^ m=^d'i3S,G; 

this is a trochoidal cuxwe, which becomes a three cusped hypocycloid for I = 2m and 
goes through the types illustrated below for different values of ^/d. For a given 
value of /3 all these ''types are present for different values of d, and are described 
about different centres owing to the comatic increase of focal length. These facts are 
well known in particular instances, and even experimentally, but as far as I can find 
they have not hitherto been expressly demonstrated generally. 

The plane at which these phenomena are found is taken at x^ = Sf\ where 
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The part depending upon d!^ represents the spherical aberration ; the part 
depending upon ^^, if it were present alone, would indicate that the images, if we 
can so call them, were found upon a sphere of curvature — K(4Gr + 4H); the 




1 = 1. m=2 






l»f, rTi«% 



1*1, m-% 



corresponding expressions for the primary focal line in place of the focal circle would 
be — K(4Gr+24H) and for the secondary focal line — K^G. These expressions are 
positive when the sphere is convex to the incident rays. 

The angular values of the radii of the focal circle and the comatic circle are 
respectively 

(1) (2) 

206265" x^4H.d/3V/' and 206265" x I^G . c^W'- 

Thus with increasing aperture {d) the focal radius increases with the first power 
and the comatic radius with the square, while with increasing breadth of field {0) the 
focal radius increases with the square and the comatic radius with the first power. 
To fix ideas we may consider the case of the parabolic reflector; here, as shown 
on p. 156, 



^,G = 0, 4G = + 1/2/', 40 = 0, 



^iH == 4 l/2j , §2^ 



OoJlJL \J, 



Hence spherical aberration and distortion are absent. S^Q = implies that the 
secondary focal line lies in the normal focal plane, while the focal circle lies upon a 
surface of curvature 1/2/^ ; and for the effects of astigmatism and coma we have the 
following table for different apertures and fields^ :™ 



Of. Poor, * Astrophysical Journal/ VII. (1898), p. 12L 
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d = //30. 



II 



15'. ' 30'. 45'. 




d = //20. 



d =//10. 



Focal radius . 0"-06 | 0"-26 i 0"-59 

Comatic radius \ 0" • 25 i 0" • 50 ; 0" • 75 
Comatic mag- ■ \ ' 

nification . 1 • 0005| 1-0005 1 -0005 
Comatic dis- 
placement . o"-5o I r'-oo r'-50 



l"-05' 0"-10 



l"-00 0"-56 



1-0005 1-0012 



2" -00 V'U 



30'. 45'. 



0" • 40 



1"'13 



0"-89 



l"-69 



60'. 15'. 



l"-57 



2" -25 



1-0012 1-00121 -0012 



2" -25 



3" -38 4" -50 



0"-19 



30'. 



0"-79 



45'. .: 60'. 



2" -25 4" -50 



1-00501-0050 



4" - 50 



9" • 00 



l"-77j 3" -14 



6" -75 9" -00 



1-0050 1-0050 



13" -50 18" -00 



The last two lines measure the same thing, the fourth representing |- ((J^H 4- 4Gr) c^^^S 
in angle, and the third l-f|-((?iH-h4Gr) (i^/3/|8/^, and owing to the relation ^Gr = ^iH 
the second line contains quantities one-half that of the fourth. Since these 
aberrations stand uncompensated, it is clear that the statement often made that the 
reflector has a very limited field is fully borne out, especially when, as is often the 
case, the ratio of semi-aperture to focal length is so great as l/lO. In this case, at 
only 30^ from the centre of the field, the light which comes from the outermost zone 
of the mirror would be spread around a little ring which was nearly a circle of 10" 
diameter, having its centre 9'' from the correct normal position for the image. 

Turn now to the other case which was proposed for discussion on p. 165, namely, 
where the original O, 0' are conjugate foci, so that H = and ^Gr = 4H. We have 
to study the delineation of any point in Oyz upon the plane O V^ or planes close 
to it. We may take the point h — const, ^ c = 0, and then make /3, y vary, so that, 
e,g,y for /3^ + y^ = const, the ray through the point (6, O) describes a cone with axis 
parallel to Ox. Let us put 

/3 = cos -v/r, y = 6 sin \fr^ 

and we have, at the plane x^ = Sf\ 

+ib [b%G -f 2be cos ^frS,G 4- e%G] -f id cos V- [b'SJI + 2bd COS V^4H + ff'S.'El 



c' + Sc' 



+ 






+iesmxlr[b%ll-{-2bdco8irSJI + e%lll{25) 



Compare these with the expressions (23) of p. 165, and we see that they run upon 
exactly the same model, but with a change of rdle in which we replace 



by 



^iGr, <^2Gr, <^3Gr, <^iH[j ^2Hj '^gH? d, <j>y j8, 

f5gH, ^sH, ^iH, Sk^G, <52G", ^iGtj 0, '\^, b. 
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♦ 

Hence it is unnecessary to work out the expressions for focal lines and the rest 
afresh since they can all be inferred without other change from what has already 
been given. 

Seidel's five conditions are usually taken as the standard form for the conditions 
of existence of a correct normal image. We may follow these and express them in 
terms of the aberrational coefficients, proceeding pari passu with the two cases :— 

0' principal focus. O, O' conjugate foci. 

(1) Absence of spherical aberration . (5iG = or ^H = 

(2) Absence of coma ^Gr = ^H = „ 4H = 4Gr = 

(3) Absence of astigmatism .... ^H = „ 4Gr = 

(4) Absence of distortion . . . . . 4H = „ ^G = 

(5) A flat field, when (2) and (3) are 

satisfied 4G = „ ^iH = . (26) 

It is of interest to consider the position occupied by the well-known conditions 
usually quoted as '' Petzval's condition for flatness of field/' and ''Abbe's sine- 
condition." 

Petzval's condition, or ^4? = 0, we see from (18) to imply 

4G--^JI = 0, 4G-4H - 0, 4K--c!iL = 0, 4K^4L = o, 

or, what is the same thing, simply 

4G = 4H and ^G = 4H 
at all distances along OV. 

If we confine attention to the two cases above, we see that in the first, where the 
emergent origin is the principal focus, G = 0, and therefore ^G = 8^B. without the 
intervention of ^^ = 0, and similarly in the second, when the emergent and original 
origins are conjugate normal foci, H = 0, and therefore ^G = 4H ; the interpretation 
of these is the same, namely, that the comatic displacement is twice the comatic 
radius — ''comatic displacement" being used to denote the expression ^cP^{8J3l-^80) 
as on p. 168 — a well-known fact, usually put in the form that in the absence of 
astigmatism the successive comatic circles have two common tangents inclined to one 
another at 60 degrees. The other term remains as the true content of Petzval's 
condition. Its interpretation may be put in different forms ; as, apart from spherical 
aberration, at the normal focal plane of any image, the longitudinal axis of the focal 
ellipse is three times its transverse axis, which is an interpretation of the expressions 
of p. 166, for 8f = 0, 

^5' = . . . cos <i>d [p^^G- + i/3^ {m + 24H)], y ^ ,.. sin <t>d i^d%G + i/3'4G], 

of the first case, or the corresponding expression of the second case ; or again, the 
distance of the focal circle beyond the normal focal plane is 2fld x the radius of the 
VOL. ccxii. — A. z 
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focal circle. This shows the connection with flatness of field, but to restrict the 

reference to curvature of the field is a misconception of the significance of Petzval's 

condition. If Petzyal^s condition holds, the result stated above is true for all normal 

image planes. 

Abbe's sine-condition for absence of coma states that if the magnification produced 

by rays passing between two conjugate foci through all zones is the same, the relation 

must hold 

sin 6^ /sin 6 = const. , 

where 6, & are the original and emergent inclinations of the ray to the axis. In our 

notation this would run 

(iS^-i/3'')/(/3-i/3') - const. ; 
or, if 

it gives 

^L/L = i(L^-l)/3^; 

now we have seen on p. 163 that the linear magnification is 

(N + M)/(L-f^L) 

and the condition this should be constant is 

^L/L-^N/N, 

and this, in accordance with p. 163, gives the conditions 

To make these agree with the sine-condition we must take & = 0, so as to remove 
^iL, S2L from the reckoning. We see then that such an assumption underlies the 
application of the sine-condition. 

I shall next show how these formulae may be applied to the numerical calculation 
of lenses. For this they are particularly appropriate if the calculations are made 
with any ordinary type of multiplying machine and not with logarithms. To the 
best of my judgment they appear to require a fraction only of the work involved in 
the equivalent complete trigonometrical calculation and, as will be shown, they are 
certainly not less accurate for telescopic object glasses. They show with remarkable 
clearness the contribution of each surface to each fault of the image. They also 
supply throughout their course a number of natural checks upon the computation 
which are searching and usually complete. 

I shall take as my example the celebrated object glass of the Fraunhofer heliometer 
at Konigsberg. This is a small lens of aperture 6*2 inches and focal length 101 inches 
which was constructed by Fraunhofer. Bessel, in describing the heliometer and 
its corrections, with his customary masterly thoroughness, calculated this lens 



MR. K A. SAMPSON: A NEW TREATMENT OF OPTICAL ABERRATIONS. I7I 



trigonometrically.'^ Later it was used as an illustration by SEiDEL,t who named his 
second condition the Fraunhofer condition, under the misapprehension that coma was 
effectively corrected in it. 

In 1889, Dr. A. Steinheil gave a particularly thorough and instructive calculation 
of its field, and assigned first, the modifications of its curves necessary to remove a 
remaining trace of spherical aberration, and next, to correct the coma.J Finally, 
FiNSTERWALDER recalculated Seidel's sums for it, using the first corrected curves of 
Steinheil. § In the following pages I shall first of all show how the calculations will 
run with my formulae, and shall return to compare them with Steinheil's results. 

The data given by Bessel are in '* lines," of which 144 = 1 Bavarian foot. To 
render the arithmetical work more compact I have increased the unit to 1000 lines, 
which brings the measures of the radii of the surfaces and the focal length of the 
whole to the neighbourhood of a unit. 

The radii and curvatures of the surfaces, with the spaces between them, are the 
following : — • 

Bo= +1193084, 

B2 = -2-996093, 

B4 = -2*936547, 

Be - - -852873, 



p^ = + -838164, 
p, == - -333768, 
p^ = - -340536, 
p^ = -1-172508, 



di — 


•006, 


<^3 = 


•000, 


d, - 


•004. 



The semi-aperture is 0'0351. 
The refractive indices he takes as 



n = m_i/mi 



•653966, 
-610083, 



l/n = 1-529130, 
l/m = 1-639121. 



We first form the normal scheme for the whole combination by writing down and 
combining the schemes that represent each surface and each space between two 
surfaces. To perform the step of combination 



9. h 



~\ 



> < 



AT, X 



> =^ < 



yg+rjk, yh-{-r]l 



> 



■ — as to which it must be remembered that light passes through {9^...} to reach 
{y ...} — we set up the number y as multipher upon the multiplying machine, multiply 
it into g and h, and place the products as above, then set up r] and multiply into k 
and Z, placing the products as shown ; this gives the top line of the combination 
completely ; then set up /c, multiply into g and h, set up X and multiply it into k and I. 

"^ * Untersuchungen,' Bd. I., p. 101. 

t * Astronomische Nachrichten,' No. 1029, p. 325. 

t K. Bayer. Akad. d. Wiss., ^ Sitzungberichte d. math.-phys. Classe,' Bd. XIX., Heft III., 1889. 

I K. Bayer. Akad, d, Wiss., ' AbJiandlungen,' Bd. XVII., Abth. III. 

z 2 
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This completes the step. The proper way to check a sequence of such combinations is 
to proceed first from left to right and then from right to left ; the final results will 
confirm one another and check the whole calculation, and the individual schemes 
arrived at will give in succession the steps across the surface (O), across the surface (o) 
and the space (l), across 0, 1 and the surface 2, and so on ; and in the reverse order 
across the surface 6, across 5 and 6, across 4, 5, and 6, and so on. All these will be 

required after. I shall indicate them with the signs 0, 01, ..., 6, 56, .... 
For the individual surfaces the values of ^ — — (l—r?.) B are 



4 =. --^(l-^n) Bo = - -412848, 
4= -^(i-^^^)B^- +1145010, 



h, := ---{l^l/n) Ba = -1 '585323, 
k^ ^ -(l"l/m)Be = - 



•545089. 



In the following arrangement the separate schemes are written in the middle 

column, that corresponding to 3 being omitted as nugatory; the calculations forwards 
are written on the left and those backwards on the right. The latter begin at the 
bottom and proceed upwards. Every figure used is recorded. 







1 • 000000 

•412848 



^ 



•653966 



1-000712 
4019 

•996693 + -006366 

• 2SS990 
• 629829 



V 



- -883819 + -997672 



06 



f 1*000000 

2^77 



01 ^ 



V, 



-> 



02 "1 



04 



997523 4- -003924 



•412848 + -653966 



1-581896 



I 



997523 + -003924 



006221 



'631298 + '999999 



2-212694 4- -993778 



< 



-997523 + -003924 

+ 1'U2174 + '00U9S 
--1' 849927 + '606287 



y 



207753 + -610780 



1 



1 • 000000 



2 



1 • 000000 



006000 



1 • 000000 



* 



1-585323 +1-529130 



4 



1 • 000000 



■?r 



+ 1-145010 + -610083 



1 



+ '006004 
+ 8781 

+ 1-000712 + -009735 

-- '001524 

+ 1 • 527096 



> 



■I 



'253990 +1-525572 

1 ' 004580 
8868 



1-000712 + -003731 

+ 1 '829225 
-- 1-588215 



U 



r 26 



253990 +1-527096 



^ 



+ -002440 



1-000000 
+ 4580 

1-004580 

- -545089 
+ 1-874814 

+ 1-329225 + -998670 



r 46 



MR. R A. SAMPSON: A NEW TREATMENT OF OPTICAL ABERRATIONS. 173 





• 997528 


•f '00892n 


5 




-> 

05 ^ 


881 


+ 2If48 


- 1 • 000000 

< 

>• 


• 004000 1 
1 • 000000 J 


•996692 


+ -006367 




~ -207753 


+ -610780 








•996692 


+ -0063671 


-> 

6 




-> 
06 ^ 


- -5^8286 

- '840582 


- '008^71 
+ 1 • 0011 Ji.2 


r 1 • 000000 

< 

L- -545089 


> 
r639121 




- •883818 


+ -997671 



1-000000 + -004000 

+ i - 6891%1 
2180 



r 56 



- -545089 +1-636941 



It may be well to repeat what these schemes imply. Take the scheme 16. If 6, ^ 
refer to any ray where it meets the tangent plane to the surface (O) after crossing 
that surface but before crossing the space 1, and h\ /3' refer to the same ray where it 
meets the tangent plane to the surface (6) after crossing that surface, then the 

scheme 16 states that 

6^ = +ro007126+ -009735^, 

/3^=:--. •253990& + 1-525572/3, 

and mutatis mutandis the same holds for c, y, c', y'. 

We conclude, from the expressions on p. 151, that for the whole combination the 
cardinal features of the normal combination are given by 

HF = -1131455 = -H'F', OF = -1128820, 6F = +1127712. 

We next work out the schemes multiplied into each of the aberrational functions o), 
as given in (ll), p. 157. The schemes [g, h, ...}, {g\ h\ ...} which respectively 
precede and follow the surface to which w refers are read at once from the computation 
of the normal system just completed. As the surfaces are supposed to be spherical, 
we have e= 1. The general arrangement is as above, and the check consists in 
forming the combination first forwards and then backwards. Again every figure is 
shown, but now the decimal places may be reduced to five. 



1-00000 



1-19308 



+ -83816 



- 1 • 00000 



i 



1 '00071 

- '01161 

+ -98910 

- '28898 
- 1 ' 82012 

-2-07410 



+ -88876 

- '00978 

+ -82903 

- 21288 
-1-52556 

-1-73844 



\ 



1 ' 00071 
- '01161 


+ '88876 
- '00978 


-98910 

- '25898 
-1' 82012 


+ -82903 

- '21288 

-1' 52556 



2-07410 



-1-73844 



1-00071 



25398 



+ -00974 
+ 1-52556 
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COo ." 



f • 99752 
I + '18780 



^ 



1-13532 

+ 2 ' 98866 
+ '41285 

+ 3-40151 



+ '00892 

- '21828 

■ -21436 

+ '01174 

~ '65897 

- -64223 



>- 



99752 



41285 



1 • 00000 



+ 2-99609 





[ 1 • 00988 
+ '18944 


+ 


'00897^ 

'22087 

•21690 


+ '003921 


1-14882 


+ -65397/ 


+ 4' 81058 
+ '59546 


+ 


'01694 

' 94828 




+ 4-90604 


__ 


•92629 




' 1 ' 00458 
+ 781 


__ 


'88580^ 

244 


~ -333771 
-1-00000 J 


1-01189 


— 


-33774 


+ 1 - 82928 
+ 2 - 99208 


— 


•44^65 

• 99866 




+ 4-32131 


, 


L-44231 



r 



4 



+ 1-14052 
+ '00880 

1-14882 

+ 1 ' 50909 
+ 8-89694 

+ 4-90603 



-« '21584 

- -00157 

- ^21691 

- -2 
~ -64187 



y 



92630 



1-00458 



+ -00244 



+ 1-32923 + -99867 



CD4 



r + - 99752 
+ -75851 

1-75103 



< 



+ 2 - 92927 
+ 2-21269 

+ 5-14196 

+ 1 ' 75108 

- 02057 

+ 1-77160 

- - 95447 
+ 8- 41708 



+ 


- 00892 

- 88848 


— 


•33451 


+ 


- 01151 

- 99879 


— 


-98228 


— 


•884S1 
• 00898 


— 


•33844 


+ 
-1 


-18284 
- 60798 



\ 



r 



[ +7-46261 -1-42559 



99752 



1 • 00000 



+ 2-93655 



1-00000 



• 54509 



+ -00392 



2-21269 + -99379 



-34054 



1 - 00000 



+ -00400 
+ 1-63694 



1 

+ 


00924 

• 76286 


1 


•77160 


< 

+ 4 
+ S 


'25180 

' 21182 


+ 1- 


46262 


+ 


00000 
'01175 


1 


•01175 

54509 
• 80696 



+ 4-26187 



+ 


• 00897 
•84240 


-— 


•33843 


+ 
-i 


• 01671 
- ff !,281 


-1 


•42560 


— - 


• 84054 
' 00400 


— 


•34454 


+ 
-1 


'18562 
• 68694 






r 



1-45132 
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O)^ .- 



^ 



99669 



20775 



+ '00636 
-h -61078 



^ 



r + '99869 
+ '24359 

+ 1-24028 

+ '85005 
+ '20775 

' +1-05780 



+ '00636 
- '71615 



- -70979 

+ '005^2 

~ '61078 

- -60536 



!. 



' + '99669 
+ '24359 


+ 


' 00636 ] 
' 71615 




1-24028 

+ '85005 

+ '20775 


__ 


-70979 

- 00542 
'61078 


r 1 - 00000 

\+ '85287 


[ +1-05780 

T- .1 


n 


• 60536 

111 ni 


■ T 1 111 



1-17251 



1 - 00000 



In the case of the first and last, the combination consisting of only two terms, the 
check calculation is a mere duplicate, and is, therefore, less searching than the others. 
The signs, in particular, should be examined to guard against a double error. 

We next form the corresponding schemes for the function i/r, again in accordance 
with the formulae (11). Owing to the occurrence of two zeroes in the scheme at the 
surface the calculation is somewhat simpler. 



\^o-- 



^ 



+ 1-19308 



+ 1 - 00000 



+ -01161 



+ '00973 



+ 1-82012 +1-52556 



+ -01161 

+ 1-82012 



+ -00973 
+ 1-52556 



1-00071 



25398 



+ -00973 
+ 1-52556 



i^2' 



^ 



i 



- 2 ' 98866 

- '41285 

-3-40151 



■X- 



- '01174 

+ '65397 

+ -64223 



r 



99752 



41285 



■X- 



-2-99609 



+ -00392 
+ -65397 



■X- 



+ 1 • 00000 



•<^ 



- '00729 

- '00101 


+ 


- 00003 
'00160 


- -00830 


+ 


-00157 


- 2 ' 98466 

- '41230 


+ 


'01173 

' 65309 


-3-39696 


+ 


•64136 


- -00731 


+: 


-00244 


-2-99208 


+' 


-99866 



00830 



+ -00157 



+ 1-00458 



+ -00244 



-3-39694 



+ -64137 



+ 1-32922 



+ -99866 
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r. 



^4^ 



99752 



2-21269 



+ '00392 
+ -99379 



^ 



•^ 



+ 5-14196 
- -02057 



- -98228 



+ -00393 



-8-41708 +1-60793 



i^e-' 



2-93655 

1 • 00000 

•54509 

•99669 
•20775 



^ 



1 



85005 

20775 



^ 



- '005^2 
+ • 61078 !^ 



85287 



+ 1 - 00000 
+ -00400 
+ 1-63694 

+ -00636 
+ -61078 

+ 1 • 00000 



' - '01172 
- '00885 


-- -02057 

•< 

^4,-79504 

~S'6220i 


-8-41708 

'■ "" -01175 
< 

--4 -80696 



-- -00005 
+ '00898 

+ -00393 

- '01884 

^1-62677 

+ 1-60793 
+ -00400 

+ 1-63694 



r 



< 



- '85005 

- '20775 



- '00542 
+ '61078 



y 



-1-05780 + -60536 



[ -1-05780 + -60536 

The schemes in the middle columns which precede and follow those belonging to 
the surface are the same for w and ^. They should be written down independently 
and read against one another to guard against errors of transcription. 

Now, for any surface (r), in accordance with (14), 

^^ = I (1 ~n,) B/ (5/ + c/), ^lr, = 1^, (/3? + y?-^/-y/) 

where 6^, c^, /3^, y^, ^\^ y\ are read from the normal schemes on pp. 172, 173 in terms of 
^0? <^oj /3o> To? the specification of the original ray. We now calculate these, noting 
that since 5^, /S^., ^\ run exclusively together, as do also c^, y^, y ^^ we need speak 
explicitly of the former only. 

We require to form from the original data |-(l— ny)B/, that is — ^^x|^B^. There 
is no check upon these values and they should be examined with care like other 
fundamental numbers. Their values are shown in the table below. 

We have then as regards m : — 



Oift 



52. 





2 
4 
6 



Co- 
efficient 

^0. 



+ 1 • 00000 
+ -99752 

ibid, 
+ -99669 



Co- 
efficient 

^0* 



Co- 
efficient 



♦ 1 1 • 00000 
+ •00392 + -99505 



ibid, 
+ • 00636 






ibid, 
•99339 



Co- 
efficient 
bol^o. 



+ •00782 

ibid, 
+ •01268 



Co- 
efficient 



•^ 



+ •00002 

ibid. 
+ • 00004 



4(1 



nr)B; 



+ • 24630 
-2-37488 
+ 1-68120 
- -23245 



(Df. 



Co- 
efficient 

b^K 



+ -24630 
-2-36312 

+ 1-67288 
- -23091 



Co- 


Co- 


efficient 


efficient 


*o/3o. 


w. 


•^ 


* 


- -01857 


- -00005 


+ -01315 


+ -00003 


- -00295 


™ -00001 



Sum 
of co- 
efficients. 



+ -24630 
-2-38174 
+ 1^68606 
~ -23387 
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The best way of forming the square, e.g., of &2? is to set up *99752 on the machine, 
multiply it into itself and into twice 392, then set up 392 and multiply it into itself 
and twice '99752, when the agreement of the middle terms is a check upon the 
operations. The last column, '' sum of coefficients in a)^," will be used below as a 
check for future work. If necessary, it may be checked by the equation, e.g. , for %, 

--2-38174 = -2*37488 X (-99752+- -00392)1 
Next, for formation of Y>r, arrange as below :— 



r. 


Pr = p r~-2- 


A-' = iSV-2^ 


• 


P r " Rr • 


J%. 


Coefficient 

^0. 


Coefficient 

+ 1-00000 
+ -65397 
+ -99379 
+ -61078 

+ -99767 


Coefficient 


Coefficient 
bo^o- 

- -53998 
-4-39790 

~ -25378 

-1-76352 


Coefficient Coefficient 


Coefficient 

boPo. 


Coefficient 

- -57232 
+ -55994 

- -61457 
+ -62230 


Sum of co- 
efficients. 



2 
4 
6 

8 


~ -41285 
-2-21269 

- -20775 

- -88382 


+ -17045 
+ 4-89600 
+ -04316 

+ -78114 


1-00000 

+ -42768 
+ -98762 
+ -37305 

+ -99535 


-17045 

4-72555 
- 4-85284 
+ -73798 


- -53998 
-3-85792 
+ 4-14412 
-1-50974 


- -94185 
+ 1-42757 
-1-32329 

- -14946 


•32698 
•76457 
-30504 
-81956 


+ -78114 


-1-76352 


- -00465 


- -98703 



The last row under /3^/— /3/ is the sum of the numbers above it and is introduced 
as a check upon the subtractions ; it is equal to /S^— i^o^. The addition of mixed 
positive and negative numbers is best done with a machine. We now have 



r. 


ij/r- 




Coefficient 


Coefficient 

bof^Q. 


Coefficient 


Sum of 
coefficients. 




2 
4 
6 


+ -05573 
+ 3*61299 
- 1 • 48032 
+ -60482 


~ -17656 
-2-94963 
+ 1-26413 
-1-23732 


- -18714 
+ -42811 

- -18747 
+ •51001 


-^ -30797 
+ 1-09147 

- -40366 

- -12249 



The multiplication by ^n^ should be checked by the help of the column '' sum of 
coefficients." It will be remembered that there is no check against setting up an 
erroneous multiplier for |-%. 

We are now ready to form §G, &c. Eef erring to (12) on p. 157 and the calculations 
above we have, for example, 

.^G= +- -98910^0+ I'14882ce)2+r77160o)4+r24028«e+ '01161^0- '00830^2 -~02057^4 + ^ 

the coefficient of 1/^2? for example, being the figure that stands in the place of G in 
the scheme formed on p. 175 for Yrg ; and similarly for SH, ^K, ^L. It is unnecessary 
to write them at length because they are shown in a more convenient place in the 
following table : — 

V OL, O CXII. ■ A. ? A A 
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We conclude from this calculation that the aberrational terms in the emergent 

ray are 

W = 5o[+-20722 (6o' + Co^)~- -00526 (Mo + Coro)-~ '00188 (A^ + yo^)] 

+ ^0 [+ '31491 (5o' + Co') + '00029(60/80 + 0070)- -00188 (^o' + yo'E 
^^' = 5o[-'21613 (V + Co^) + -37086 (60/80 + ^oro)- '^5647 Oo^ + yo')] 

+A[--'09612(6o^+CoO--87730(6o^o + c,yo)-~-0036lOo' + yo^)] . (27), 

with similar expressions for 8c\ Sy if we replace 60, ^80 outside the brackets by Co, yo- 
We may without loss of generality put yo = 0, and in what follows this shall 

be done. 

The twelve coefficients above are not independent. We have seen that they must 

satisfy seven relations, and we shall now verify that they do so. 
We have, from p. 161, 

5|j = ^0 + ^^2 + ^4 + ^^6- 

h 

- -41285 
-1-58532 
+ 1-14501 

- -54509 



and thus 



53K 



SiL 



- -00526- -62982 

- -00376- -00029 
+ -37086 4- -19224 
-1-51294 4- -87730 



Term in % 






- -41285 






-- 1-03675 






+ 1-14501 






~~ -33255 


= * 




~- -63714 = 




- -63508 


^a = 


"-63503 


- -00405 


$H = 


- -00406 


+ -56310 


glK = 


+ -56312 


- -63564 


qjL = 


- -63566 



Next for the relations (20) ; we have, noting that N = 1, 



s. 


+ LS,G. 


- K5,H. 


~ HS,K. 


+ G5,L. 


- 4N. 


1 


+ -41347 


+ -55665 


+ -00275 


- -19160 


+ -78127 


2 


- -00525 


+ -00026 


- -00236 


- -87440 


- -88175 


3 


- -00376 


- -00332 


+ -00963 


-- -00720 


- • 00465 



NK2 = + -78113 



NKL - 



88176 



N(L2-1) = --00465 



The comparisons seem to point to a small accumulated error in ^^H. 
If we transfer the origin for emergent rays to the focal plane, x' = 4-1*12771, we 
have finally 

a6' =: 6 [--•03651(12 4- -41296 cl^ cos 9^- -85496^^; 

+ ^[+ -20651 #- -98905^^0.08 ^--00595/3^] 

Sc' -= €[--03651c?2 4--41296d:/3cosf--85496^'] 

2 A 2 



# 



(28), 
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where the subscript (O) is dropped,// is taken as zero, and b ^ d cos (j>, c ~ d sin (p^ 
as on p. 165. 

These give the aberrations of the lens at its principal focus. 

The corresponding normal scheme is 

V = ^ +1'131453^, 

- /3' := ^•883818&+ '997671/3. (29) 

In order to fix ideas, compare (28) with the case of a parabolic reflector of the 
same focal length, given on p. 167, for which we get 

8V =:h\_ ^ +-44191&/3+ '^ ] + ^[ + '22095c?2-l-00000fe^+ ^ ]; 

we see that there is a close resemblance, except for the value of ^G, so that the two 

hardly differ in any sensible way, except in the curvatures of the fields. It will then 

cause some surprise that Seidel concluded that the Fraunhofer glass was free from 

coma which is so marked in the reflector. It was, in fact, a misapprehension, as the 

diagrams given by Steinheil and by FiiSfSTERW alder sufficiently demonstrate. 

Seidel's argument presents an interesting feature.! He puts together the four 

components of his sum S (2), 

+ 0,412 

-12,672 

+ 13,454 

— 1,662 



S(2) = - 0,468 

and draws his conclusion from the approximate balance, within one-thirtieth, of the 
large positive and negative members. It is evident, however, that this amounts to no 
more than saying that the two internal surfaces nearly annul one another. But the 
point I wish to make is that these numbers are in fact the same as those found 
on p. 178 above. If we transfer to the principal focus by adding /'^K to ^G, we have 
for 4Gr 

(0. ^. S2O, S2G x/. 

* - -36445 = - '36445 + -4125 

-12406 +11-32385 +11-19979 -12-6721 

+ -13396 -12-02515 -11-89119 +13-4543 

-•00718 + 1-47599 + 1-46881 - 1-6619 



+ -41296 + 0-4672 

The connection does not appear to be so close in the case of others of Seidel's sums, 
but it is interesting to notice this common ground. 

Let us now compare my calculations with those of Steinheil. First as to 

t ' Ast. Nach.,' No. 1029, 326. 
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spherical aberration. Take rays parallel to the axis, that is, /3 = 0, and consider the 
focus where they unite for impact upon the original plane at fractions 0, l/3, 2/3, 
and 1 of the semiaperture, i.e., for 

d = O'OOOOO, 0^01170, 0^02340, O'OSSIO. 

For any of these the distance of the point from the last surface (6) is 



or 
and 



G/K[l+idHS,G/G-S,K/K)l . , (30) 



G = +^996692, K = -•883818, ^.^iG = +"20722, ^S^K 

:. |-[^iG/G-^iK/K] = +"20791-^24454 = --03663. 
Hence the rays meet at the following points along the axis : — 



•21613. 







Steinheh., p. 417. 


Axial 


. 1127'-712 


1127'^712 


1/3 semiaperture . 


•706 


•706 


2/3 


•689 


•687 


1 


•662 


•659. 



. . (31) 

In these and the following comparisons the unit of length has been brought back 
to 1 Hne by multiplying by 1000, to preserve Steinheil's numbers unchanged. 

In consequence of this residue of spherical aberration the best setting for focus at 
the middle of the field is not the axial focus but a point within it. Steinheil takes 
this point at 1127*670, following presumably the theory of Bbssel, which gives a 
position for the greatest apparent concentration of light that is slightly within the 
least circle of aberration (1127*672).^ 

Adopting the corresponding point, which allows for the slightly smaller aberration 

shown by my numbers, and multiplying by 10"^ to bring the units into agreement 

with formula (28), we see, in accordance with p. 165, we must include with W of 

p. 179 the term 

+ K5/'.6= --8838 X "^'0000400x6 = +*000035356, 

with a corresponding term for 8c^ in terms of c. 

The diameters of the image-disc in the focal plane and at this setting are 
respectively the corresponding extreme values of 8h\ doubled, or 



Steinheil. 



0^-00316 
0^-00067 



0^-00071. 



■^ Bessel, he. cii, p. 104. 
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We take next the oblique rays in a plane through the axis, that is, we take c = 0. 
The rays considered are taken at an angle of 48^ with the axis, following Bessel 
and Steinheil; hence /3 = tan 48^ = '01396353, and we take in succession 
6 = cZ= +'03510, +'02340, +*01170, -'01170, -'02340, -'03510. The final 
numbers given below have been multiplied by 1000 in order to compare with 
Steinheil. 

The central ray meets the chosen plane at a distance from the axis 

(H + L^/0/3 = ( + l'131453-'000040) x j3 = '01579852. 

Giving as before the calculations in full, the formulae (28), supplemented by the 
term K.Sf^ . &, give the following : — ■ 



The comma is placed between the 7*^ and 8"' decimals. 
Hence we have the following, with unit 1 line :— 



b. 


8b\ 


h' + Sb\ 


Steinheil, p. 419. 


+ 35^-1 


+ 91 323 ~ 


: - -00232 


15^-79628 


15^-79622 


23 -4 


~ 38- 295 


333 


•79519 


•79528 


-f 11 -7 


- 81- 188 


269 


•79583 


•79587 


-0 


0- 2 


2 


•79850 


•79852 


-11 -7 


+ 238+ 263 


+ 501 


•80353 


•80349 


23 -4 


+ 669+ 608 


+ 1277 


•81129 


•81130 


-35-1 


+ 1329 + 1031 


+ 2360 


•82216 


•82212 



b. 


+ KS/'. 


+ i3^Gd^. 


+ S^Gb/3. 


+imf^^ 


Coeff. b. 


+ ^8,RdK 


+ mhl3. 


+ h^-MI3^ 


Coeff. ft. 
-2314,7 


+ •03510 


+ 353,5 


~~ 449,8 


+ 2023,9 


- 1667,2 


+ 260,4 


+ 2544,2 


- 4847,3 


-11,6 


•02340 


ibid. 


-199,9 


+ 1349,3 


ibid. 


- 164,3 


+ 1130,8 


- 3231,5 


ibid. 


-2112,3 


+ •oino 


ibid. 


" 50,0 


+ 674,6 


ibid. 


~ 689,1 


+ 282,7 


- 1615,8 


ibid. 


- 1344,7 


- -01170 


ibid. 


~ 50,0 


- 674,6 


ibid. 


-2038,3 


+ 282,7 


+ 1615,8 


ibid. 


+ 1886,9 


•02340 


ibid. 


-199,9 


^ 1349,3 


ibid. 


-2862,9 


+ 1130,8 


+ 3231,5 


ibid. 


+ 4350,7 


- -03510 


ibid. 


-449,8 


- 2023,9 


ibid. 


- 3787,4 


+ 2544,2 


+ 4847,3 


ibid. 


+ 7379,9 



(32) 



There is a slight discrepancy for 6 — +23 '4. 

Steinheil now considers the rays which do not meet the axis. Fig. 6 is 
taken from his memoir and shows the object-glass on reduced linear scale. He 
divides the object-glass into three rings, and computes all the rays which impinge 
upon it at an angle of 48^ with the axis, at the points indicated in the figure. The 
rays 2, 10, 18, 1, 22, 14, 6 are those just given ; of the remainder, those upon the left 
may be written down from symmetry from those upon the right, so that he computes 
in all nine independent rays which do not meet the axis. 
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We derive these as follows. We have throughout ^ unchanged : — 

For the rays (4), (12), (20), 6 = 0, c == d; 

,, . ,, (3), (11), (19), 6 = c = c?sin45° = d-x70711, 
„ (5), (13), (21), -b =:o = c?sin45° = c?x7071L 




Fig. 6. 

Hence the additional calculations required run as follows. 

In the calculation of the coefficients replace the columns 4Gr6/3, d^h^ by the 
following : — 




In the rays (4), (12), (20) replace these by zero. 
Therefore 



rv— ' -J„ J^JL JL^ -1 — I. JL-f JL JLwl • 



Ray. 


Coefficient 
h ox c. 


Coefficient 




- m\ 


8c\ 


3 


- 332,4 


- 895,0 


- 82-125 = 


-^•00207 


-^•00082 


11 


- 559,5 


-1165,8 


- 93-163 = 


- -00256 


- -00093 


19 


- 886,7 


™ 871,4 


- 73-122 = 


- -00195 


- -00073 


21 


- 1840,7 


+ 1413,6 


+ 152 + 197 - 


+ -00349 


- -00152 


13 


-2467,7 


+ 3404,2 


+ 408 + 475 - 


+ -00883 


- -00408 


5 


--3194,6 


+ 5960,2 


+ 793 + 832 = 


+ -01625 


~ -00793 


4 


- 1763,5 


+ 2532,6 


+ * +354 = 


+ -00354 


- -00618 


12 


-1513,6 


+ 1119,2 


+ ^ +156 = 


+ -00156 


- -00354 


20 


- 1363,7 


+ 271,1 


+ * +038 = 


+ -00038 


- -00160 





W. 


-^-00204 






•00251 
•00191 




+ 


•00351 




+ 


'00879 




+ 


•01605 


. 


+ ' 


•00357 




+ • 


00160 




+ • 


00041 



8c'. 



-^•00080 

- -00092 

- -00070 

- -00151 

- -00410 

- -00758 

- -00627 

- -00353 

- -00162 



(33) 



There is a slight discrepancy in the ray (5). 
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In considering what discrepancies may be expected, we have to recall that the 

method developed in the preceding pages omits terms of the fifth order, which may 

amount to, say, 

coe^xdx -000001. 

Taking d ™ 35^'l, for unity as coefficient, we should have an error of 4 units in 
the last place retained above. We have no means of saying what the coefficient may 
be, but it is clear that it may affect the last digit. Yet I believe that these 
calculations are not only very much easier, but also more correct than the trigono- 
metrical ones, for though the formulae for the latter are exact, the number of 
operations they require is very large. Thus, for each ray which meets the axis, there 
are fully 50 operations of which at least one-half consist in taking out a logarithm 
or an antilogarithm with seven decimal places ; for each ray which does not meet the 
axis the work is rather more than four times as great, 

Steinheil has calculated seven of the former rays and nine of the latter. 

The controls that exist are of the most meagre description and give little help in 
locating an error. But, even if the whole is done in strictest accordance with the 
tables, at any step an error may be introduced which falls only short of half* a unit 
in the last place. Thus, in the rays which do not meet the axis, an irremovable 
accumulated error of 10 or more units could cause no surprise, and for this reason the 
trigonometrical method loses any advantage over the formulae given above which it 
might claim from resting upon exact formulae. The differences under discussion are, 
however, minimal, since 550 units in the last decimal place only amount to 1 second 
of arc. 

But pursuing the question a little further I believe, in spite of the evident care 
with which the whole of Steinheil's calculations have been carried through, that the 
comparison above shows that a small error has crept in in respect to ray (5). 

If we take the general agreement as showing that the trigonometrical calculation 
does in fact bring in no terms of the aberrations beyond the 3rd order, we can readily 
analyse Steinheil's numbers in more than one way so as to derive the coefficients 
(?iG, . . . from them. 

Take the formulae (24); on the outer ring = 0, 45°, 90°, 135°, 180°, correspond 
respectively to the rays 2, 3, 4, 5, 6 ; thus we have 

d^/34G = ^c'3- ^c\ = i {Sb', + 56'e) - ^^'4 

KdSf+^d%G+^dl3' (4G + 24H) - = -l^{Sb^,^§V,) - USV,--^b^,) 

d'^ {SJl + S,G) + ^%R - Sh', + Sb\ - i {SV, + §b',) + W, 

P 03X1. = 06 1 . . » » . r • , . . , . i"4j 
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From these we get at once 



and thence 



d^^4G, ^^4H, and d^%II, 
KdSf+ ^d\G + id/3%G and d'^SJl. 



Also for the case /3 = 0, 



Sb^, = KdSf'-^id%G ; 



thus we get d/3^SgG, and when the adopted value of Sf^ is used, d\G also, which 
completes the solution. 

We see that we can use the rays (A) — 3, 5 exclusively, or (B)— 2, 4, 6 exclusively. 
Making separate determinations by these roads, 









(A.) (B.) 




d'/3S,G 






+ •00678 +^00709 




KdSf 


+ ld%G+ldl3%G . 




-•00593 -^00627 




KdSf 


+ ld%G+id^^SsG+ 


24H) . . . 


-•01279 -•01295 




d'/3 {§,11 +S,G)+/3%Il. . 




+ •01401 +^01422 






p o^tL .... 


• • • • 


. . •ooooo 






2Kd§f'+d%G . 


• • • • 


, . -•00071 






2KdSf' . . . 




+ •00260 




Hence 












(A.) (B.) 




(A.) (B.) 


p. 179. 


d\G . . 


-•00331 


oiG . 


-•0766 


■0730 


d'/3S,Q . . 


+ •00678 +^00709 


02^ • 


. + ^3941 + ^4121 + 


•4130 


dl3%G . . 


-•01115 --01183 


^sG . 


1^6292 1^7285 


r7 09 9 


d'issja. . . 


+ •00723 + -00713 


(^jH . 


. + ^4202 + ^4144 + 


•4130 


d^'S^a . . 


-•00686- -•ooees 


4H . 


. -r0023 - ^9760 - 


•9891 


^%R . . 


•00000 


4H. 


•0000 - 


•0060 



There is no doubt, from the checks on p. 179, that the numbers put in the last 
column for comparison are correct to the last digit, and we see that the numbers (A) 
which rest upon the ray (5) are decidedly less consistent with them than the 
numbers (B) which do not. 
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